Surface and Image-Potential States on the MgB 2 (0001) Surfaces 
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We present a self-consistent pseudopotential calculation of surface and image-potential states on 
MgP>2(0001) for both B-terminated (B — t) and Mg-terminated (Mg — t) surfaces. We find a variety 
of very clear surface and subsurface states as well as resonance image-potential states n=l,2 on both 
surfaces. The surface layer DOS at Ef is increased by 55% at B — t and by 90% at the Mg — t 
surface compared to DOS in the corresponding bulk layers. 
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The discovery of the superconductivity in a simple 
metal polycrystalline compound MgB 2 with the critical 
temperature T c ~ 39 K Q has generated an explosion 
of research activity in studying the mechanism of the su- 
perconductivity and properties of this compound [p|-p5| . 
For instance, the superconductivity gap has been mea- 
sured by both bulk sensitive methods and surface 
sensitive techniques. Compared to bulk measurements 
the surface sensitive experiments, namely scanning tun- 
neling spectroscopy (STS) and point-contact exper- 
iments fllcl , give generally a smaller energy gap varying 
in the surface region |Q,|l0|. It may be caused by two 
effects: surface contaminations or/and desorder and by 
changing the electronic structure at the surface. Qualita- 
tively different STS spectra obtained by different groups 
on polycrystalline MgB 2 pellets and films reflect different 
surface contaminations and microstructure of the sample 
surfaces. However, for a single crystal a possible change 
of a high density of states (DOS) at the Fermi level, Ep, 
high phonon frequencies and strong electron-phonon in- 
teractions can also lead to a change of the energy gap 
and T c at the surface. Very recently two groups have 
announced obtaining single crystalline MgB 2 with edge 
angles of 120 degrees |26|,^7]]. These works open up new 
prospects for experimental studies of the surface super- 
conductivity in MgB 2 . 

Due to strong covalent interactions within B planes 
0,^2|j22[|2|l the (0001) termination oiMgB 2 is supposed 
to be more favorable. However nothing is known about 
atoms which form the topmost layer of MgB 2 (0001). The 
study of the (0001) termination of other metal diborides 
which also have crystal structures of the AIB 2 type have 
shown that some metal diborides (TiB 2 , HfB^) are ter- 
minated by metal atoms |^8|,^9| while the topmost layer 
of TaB 2 is formed by the graphitic boron layer |3(| . Here 
we report on ab initio calculations of the electronic struc- 
ture of the MgB 2 (0001) surface for both types of termi- 
nation. In order to assess the effect of the surface re- 
laxation on surface states we have computed the surface 
electronic structure for the ideally bulk terminated crys- 
tal as well as for contracted and expanded by 6% the first 



interlayer spacing. 

The bulk electronic structure of MgB 2 Jl|@||,§| 
leads to an unconventional bulk states projection with 
very wide absolute and symmetry energy gaps (Fig. 1) 
which support a variety of surface states and give an addi- 
tional contribution to crystal reflectivity at energy inter- 
val just below the vacuum level where resonance image- 
potential states arise. The surface and image states are 
of crucial importance for the description of the surface 
dynamical screening, electron (hole) excitations, and su- 
perconductivity at MgB 2 surfaces. We show that for the 
Mg-terminated (Mg — t) surface the surface states contri- 
bution nearly doubles the surface DOS at Ep compared 
to the bulk Mg layer DOS. For the 5-terminated (B — t) 
surface the surface state contribution increases the sur- 
face DOS at Ep by 55% compared to that in bulk. Spe- 
cial attention in the paper is focused on image potential 
states. We find that Mg- and B-layers possess distinct 
reflectivity that leads to different localization of image 
state wave function in bulk region. 

The calculations of charge density have been performed 
within the self-consistent local density-functional plane- 
wave method by using a supercell of 13 atomic layers and 
7 layers of vacuum [|lj. This supercell is big enough to 
ensure a good description of both surface and bulk states. 
Experimental values of lattice constants a= 5.8317 a.u. 
and c= 6.6216 a.u. used in the evaluation have been 
taken from Ref. 0]. The 13 layer slab representing the 
B — t (Mg — t) surface consists of 7 B(Mg) layers alter- 
nating with 6 Mg(B) layers. 

As the LDA potential does not describe the correct 
asymptotic potential behavior in the vacuum region we 
modify it by retaining the self-consistent LDA one for 
z < Zi m , where Zi m is the image plane position, and re- 
placing it in the vacuum region for z > Zi m by V(z) — 
{exp[-X(z - Zi m )] - l}/[4(^ - Zim)]- The damping pa- 
rameter A is a function of (x, y) and is fixed by the re- 
quirement of continuity of the potential at z = Zi m for 
each couple of (x,y). With the use of the self-consistent 
charge density obtained for a 13-layer slab we have con- 
structed charge density for a 31-layer slab inserting 18 
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bulk layers into the center of the slab. Vacuum space 
was increased from 7 layers to 21 ones. This vacuum 
interval is enough to accurately describe the n=l and 2 
image states. Finally the LDA potential was generated 
for this new supercell with a correct image tail in vacuum 
space. 

In Figs. la and lb we show the calculated projection of 
bulk band structure onto the surface Brillouin zone to- 
gether with surface states for B — t and Mg — t surfaces, 
respectively. The light grey areas show the n projected 
bulk states and the grey ones indicate the a states. A 
remarkable feature of the bulk states projection is the 
presence of two wide absolute energy gaps. The lower 
gap separates the s-bulk bands and p z bands of boron, 
the upper gap crossed by the p XiV bulk bands of B is lo- 
cated in the vicinity of the Fermi level, Ep. The B—t sur- 
face (Fig. la) has 4 surface states strongly localized in the 
topmost boron layer and 3 subsurface states. All these 
states show energy dispersion which repeats that of bulk 
bands. Two surface states degenerated at the T point 
are of the p x ^ y symmetry {a states). They split off from 
bulk states of the same symmetry by 0.45 eV and have 
energy of 1.23 eV relative to Ep. Their charge density 
is completely localized in the topmost layer (Fig. 2). One 
can consider these states as two-dimensional quantum- 
well states due to their extremely strong localization in 
the z direction: they decay into the bulk much faster than 
do conventional surface states which are characterized by 
a smooth exponential decay. Another surface state with 
energy of -2.74 eV is of the p z symmetry (n state), 75% of 
the state is concentrated in the three surface atomic lay- 
ers and in the vacuum region (Fig. 3). The lower surface 
state is of s symmetry and splits off from bulk states by 
0.4 eV, 70% of the state is localized in the topmost layer. 
The subsurface states with energy of 0.35 eV degenerated 
at T are localized in a few subsurface B layers with 40% 
of the state concentrating in the second B layer. The 
third subsurface state is located at the bottom of the s 
bulk boron states. 

The Mg — t surface shows distinct electronic structure 
compared to the B — t one. In particular, the Mg occu- 
pied surface state of the s — p z symmetry with energy of 
-1.94 eV appears at the V point. Its charge density dis- 
tribution localized mostly (65%) in the Mg surface layer 
and in the vacuum region is shown in Fig. 4. The ori- 
gin of this state can be understood from a simple charge 
transfer picture. In bulk the Mg atom donates two va- 
lence electrons to the adjacent B planes thus moving up 
to E > Ep all the Mg bands. In the surface layer the 
Mg atom donates one electron to the subsurface B plane 
while another electron forms the dangling bond (s — p z ) 
occupied surface state. 

Unoccupied subsurface states with energy of 0.36 eV 
degenerated at T are formed by the subsurface B layer, 
70% of the state is concentrated in the layer. At energy 
~ -12.3 eV there also exists a subsurface resonance state 



generated by the B layers. 

In Fig. 5 we show the calculated surface layer DOS for 
both the B—t and Mg—t surfaces and compare them with 
the corresponding central layer DOS. In the B — t surface 
the surface DOS at Ep which also includes the vacuum 
region is by 55% higher than the central B layer DOS. In 
the Mg — t surface the surface DOS at Ep is higher than 
the central Mg layer DOS by factor of 2. Both these 
results favour the higher surface critical temperature T c s 
compared to that in bulk. 

It is less known about phonons on the Mg£>2(0001) 
surface. There normally exist surface phonon modes on 
metal surfaces with slightly smaller frequences compared 
to those in bulk S). In bulk MgB 2 the in-plane boron 
mode Ei g is responsible for strong electron-phonon in- 
teraction lyim. Because of its in-plane character one 
can expect that the vibrational frequences and atomic 
displacements of this mode in the surface or subsurface 
boron layer will be similar to those in bulk. Therefore 
one can expect very similar or even higher T c s compared 
to T c . 

Image states fall in a group of surface states which 
are linked to the vacuum level and located relatively far 
from the surface. The calculated work function which 
fixes the vacuum level relative to Ep was obtained to be 
of 6.1 eV for the B — t surface and 4.2 for the Mg — t one. 
Similar to simple and noble metal surfaces ]33| the wave 
function maximum of the n=l image state on MgB2 is 
located at ^6 a.u. beyond the surface atomic layer for 
both surfaces. In Figs. la and lb we show the calculated 
n=l,2 resonance image states. Nothing is known about 
the image plane position on MgB 2 and we varied Zi m for 
both terminations within the 2.0-3.5 a.u. interval beyond 
the surface layer. This variation leads to E\ = — 0.9±0.15 
eV and E 2 = -0.25 ± 0.05 eV for the B — t surface as 
well as to Ei = -1.1 ±0.15 eV and E 2 = -0.30 ± 0.05eV 
for the Mg — t surface, the error bar includes the energy 
dependence on the Zi m position. The energies obtained 
are rather similar to those for the n=l,2 resonance image 
states on simple metal surfaces . 

The resonance image states are mostly degenerated 
with magnesium bulk states. The interaction between 
the image states and the Mg bulk states results in a dif- 
ferent reflectivity of B and Mg layers and a different 
behavior of the image state wave functions in bulk. The 
amplitude of these wave functions is significantly larger 
in magnesium layers than in boron ones. This behavior 
of image states is specific for MgB 2 due to its peculiar 
bulk electronic structure and was not found for simple 
and noble metals (3^] . 

It is known that the relaxation of closed-packed sim- 
ple metal surfaces J34| is relatively small: the contrac- 
tion/expansion of the first interlayer spacing is < 6%. 
We have inspected the dependence of the surface elec- 
tronic structure by computing slabs with the contracted 
and expanded first interlayer spacing by 6% for both ter- 
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minations. We have found that these relaxations lead to 
change of the surface state energies within 0.1 eV and 
to small change of the surface DOS at Ep. The change 
of the n=l,2 image state energies is significantly smaller 
than the error bar. 

In conclusion, we have performed self-consistent pseu- 
dopotential calculations of the surface electronic struc- 
ture for the B — t and Mg — t surfaces of MgB^. We have 
found a variety of surface and subsurface states as well as 
two resonace image states on both surfaces including an 
unoccupied quantum well state of the p x . y symmetry on 
the B — t surface. Due to very clear surface character of 
these states the Mgi?2(0001) surfaces provide good opor- 
tunity to test the theoretical results by measuring the sur- 
face electronic structure by different spectroscopies such 
as photoemission, inverse photoemission, time-resolved 
two-photon photoemission, and scanning tunneling spec- 
troscopy. The obtained higher surface layer DOS at Ep 
favours the higher critical temperature compared to that 
in bulk. This is inconsistent with recent STS experiments 
which have shown an opposite trend (| We attribute 
this discrepancy to contaminations and desorder on poly- 
crystal sample surfaces. 
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FIG. 1. Projected bulk band structure of MgB A together 
with the surface and image-potential states for B-terminated 
(a) and Mg-terminated (b) surfaces. The light (dark) grey 
areas represent the -k (<j) projected bulk states. Dotted areas 
indicate the magnesium projected states. Thick solid lines 
depict surface state which are mostly localized in the topmost 
B layer. 



FIG. 2. Charge density distribution of the unoccupied 
boron surface (quantum well) state at F in the (1010) plane 
for the B-terminated surface. Small (large) filled circles indi- 
cate the B(Mg) atom positions. 



FIG. 3. Charge density distribution of the occupied boron 
p z surface state at V in the (1010) plane for the S-terminated 
surface. 

FIG. 4. Charge density distribution of the occupied mag- 
nesium p z surface state at F in the (1010) plane for the 



Mg-terminated surface. 



FIG. 5. The calculated surface layer DOS at Ef (solid 
lines) for both B- and Mg-terminated surfaces. Dashed lines 
show the difference between the surface layer DOS and the 
corresponding central layer DOS. 
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